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We report studies of photon-stimulated desorption (PSD), also known as light-induced atomic desorption 
(LIAD), of sodium atoms from a vacuum cell glass surface used for loading a magneto-optical trap (MOT). 
Fluorescence detection was used to record the trapped atom number and the desorption rate. We observed a 
steep wavelength dependence of the desorption process above 2.6 eV photon energy, a result significant for 
estimations of sodium vapor density in the lunar atmosphere. Our data fit well to a simple model for the loading 
of the MOT dependent only on the sodium desorption rate and residual gas density. Up to 3.7 x 10 7 Na atoms 
were confined under ultra-high vacuum conditions, creating promising loading conditions for a vapor cell based 
atomic Bose-Einstein condensate of sodium. 

PACS numbers: 34.50.Cx, 37.10.De, 37.10.Gh, 37.10.Vz, 68.43.Tj, 95.30.Dr, 96.12Jt, 96.12.Ma 



L INTRODUCTION 



LIAD (light-induced atomic desorption) is a new tool for 
producing laser cooled samples of K, Rb and Cs JH-Ht]. 
Its primary appeal is the non-thermal nature of the des- 
orption process J^,01, which allows for fast switching- 
the gaseous vapor can be immediately quenched upon 
shut-off of the light. Even in tiny quantities, residual va- 
por can dramatically reduce the efficiency of evaporative 
cooling to Bose-Einstein condensation (BEC), and LIAD 
has been shown to mitigate this problem [8]. Other so- 
lutions to the background vapor problem involve loading 
from differentially pumped secondary sources such as a 
dual MOT or slow atomic beam. Typically these are more 
cumbersome approaches in comparison to a single vapor 
cell. 

In this paper we report the use of LIAD to load a sodium 
MOT with up to 3.7 x 10 7 atoms in a pyrex vacuum cell. 
Moreover, we demonstrate that such a MOT is compat- 
ible with vacuum limited lifetimes in the several second 
range. We present measurements of the wavelength de- 
pendence of the LIAD effect, showing that it increases 
dramatically in the blue to ultraviolet region of the spec- 
trum, consistent with a photoelectric threshold. 

Photon-stimulated desorption (PSD) studies of neutral 
sodium are particularly important for a better understand- 
ing of the lunar surface and atmospheric composition. 
Sunlight irradiating the silicates present on the moon's 
surface creates a significant sodium and potassium vapor 
in the tenuous atmosphere, causing the moon to glow due 
to emission from the same vapor Optically induced 
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FIG. 1: LIAD assisted loading of a sodium MOT. Atoms were 
captured from a transient vapor created by LIAD. Shown are 
MOT loading curves for different light intensities using a blue 
high-flux LED (455 nm): 0, 0.9, 2.6, and 4.0mW/cm 2 . Simple 
exponential loading curves fit each data set and the inverse time 
constant y extracted from the fits is presented. We observe that 
the number of trapped atoms in the steady state increases mono- 
tonically with the LIAD intensity, saturating at high intensities. 
In the inset graph, the MOT loading rate is plotted as a function 
of the LIAD intensity, presenting a linear dependence, in very 
good agreement with a simple model described in the text. 



desorption studies have been performed on small sodium 
clusters lloll . Sapphire lfl"j|and SiOjjO surfaces, as 
well as thin organic films yl 0, [T3l4l5ll . Our measure- 
ments on atomic sodium are the first to use a MOT to 
probe PSD, and they may open up the possibility of us- 
ing laser cooling techniques for chemical identification 
of lunar samples. 

Apart from extraterrestrial applications, we believe our 
results will contribute to the widening use of vapor cell 
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technology for producing better controlled Na MOTs in 
the range of tens of millions of atoms Bio, Il7ll . Indeed, 
a recent preprint reports a vapor cell based sodium BEC 
by rapid evaporation in a dipole trap, demonstrating that 
Zeeman slower technology is not required |18] . The 
sodium atom has many favorable properties for Bose- 
Einstein condensation, including a large ratio of elastic to 
inelastic scattering rates 11911 . Its light mass contributes 
to large quantum mechanical effects, including the Bose- 
Einstein phase transition temperature, the separation of 
matter wave interference maxima, Josephson tunneling 
probability and atomic recoil shift [20]. These phenom- 
ena are all of fundamental importance in addition to their 
significance for building practical atom optical devices. 



II. EXPERIMENT 



We used light from a single mode dye laser (Coherent 
899, locked to near the F — 2 — > F' = 3 hyperfine 
components of the D2 line of sodium at 589 nm by sat- 
uration spectroscopy) to make a standard retro-reflected 
three beam magneto-optical trap B2lll . Prior to entry into 
the vacuum chamber the trap laser beam passes through 
a 1712 MHz electro-optic modulator (New Focus model 
4461), which introduces about 20% sidebands on the 
F — 1 — > F' = 2 repumping transition. Typically, 
250 mW of total laser power was used from trapping, 
split into three independent beams of 20 mm diameter 
(laser power and beam diameters chosen to optimize the 
MOT), and traveling along the X,Y and Z directions. A 
pair of quadrupole coils created the magnetic field gradi- 
ent of 20 gauss/cm (at I = 15A) along the vertical direc- 
tion. 

The vacuum chamber comprises a Pyrex cell of dimen- 
sions 45 x 45 x 150 mm attached to one face of a stain- 
less steel 2 3/4" ConFlat cube (MDC). The cube is also 
attached to a 20 1/s ion pump and titanium sublimation 
pump, and could be evacuated through an all-metal seal 
valve. The titanium sublimation section was added to 
achieve UHV conditions (pressure below the lower limit 
of the vacuum gauge, < 1 x 10" 11 torr) and the experi- 
ments with HV (pressure as 6x 10~ 9 torr) were done with- 
out this additional pump. The overall dimensions of the 
chamber were 50 x 55 x 52 cm, rendering it light and 
easily removed for modifications or maintenance. 

The sodium source was a set of alkali dispensers 
from SAES (SAES Alkamax metal dispensers, model 
NA/NF/1.5/12 FT 10+10, active length 12 mm, thresh- 
old current 6 A). Copper wires extending from a vacuum 
feedthrough were used to run current to the dispensers 



as well as to secure them about 1 cm inside the cell. 
The actual distance between the MOT and the dispenser 
could be varied from 15 mm to 70 mm by simply trans- 
lating the vacuum chamber relative to the MOT center. 
No significant variation of the number of trapped atoms 
was observed, and therefore during the experiments this 
distance was fixed around 40 mm. For the HV experi- 
ments, one dispenser was fired at 5 to 7 A during exper- 
imental runs to fill the cell with sodium. The MOT coil 
current was used to turn the trap on and off. For UHV, 
we first fired the same dispenser at no more than 3 A for 
30 minutes to coat the cell walls and waited for a few 
hours until the pressure recovered. In this case, a load- 
ing sequence was initiated by turning on the LIAD light 
a few seconds before the coil current was on in order to 
allow the LED to thermally equilibrate. In both cases, 
fluorescence from the MOT was recorded by a calibrated 
photodetector whose output current is proportional to the 
atom number. The photocurrent was recorded on an os- 
cilloscope and comprised all of the data recorded and pre- 
sented here. The calibration between photodetector cur- 
rent I p h and atom number N can be written as 12111 



N = 



F sc he/ 'A D. 



(1) 



where r sc (7) = 



s + 4(5 2 /r 2 ) is the laser-intensity 
dependent scattering rate, s = I/I sat is the satura- 
tion parameter with I sal = 6.3mW/cm 2 for Na. 6 = 
-15 MHz was the detuning and F = 9.8 MHz is the nat- 
ural linewidth. The other parameters are the responsivity 
of the Si photodector9i, = 0.38 A/W atA = 589 nm, and 
Q = 0.1 sr the solid angle covered by the photodetector. 

With the exception of Figure|5] all data was taken by irra- 
diating the cell walls more or less uniformly using a com- 
mercial blue LED (Optotech OT16-5100-RB, A LIAD = 
455 ± 20 nm, and 480 mW measured maximum output 
power). 



III. DATA ANALYSIS 

A typical set of loading curves is shown in Figure Q] and 
demonstrates the photon stimulated desorption of sodium 
from the cell walls. At t = the MOT current is turned 
on. Subsequently, the number of trapped atoms increases 
and reaches a steady state whose value depends mono- 
tonically upon the LIAD intensity. These curves can be 
analyzed in terms of a simple rate equation model for the 
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number of trapped sodium atoms N(t) [2 
N(t) =R- (y Na + y bg )N(t) 



(2) 



where R is the MOT loading rate, -y-^a and yb e are the 
trap loss rates due to two-body collisions between the 
captured atoms with fast untrapped Na atoms and with 
background gases, respectively. We have neglected intra- 
trap light-assisted collisions which appear as an addi- 
tional quadratic term in the sodium density and did not 
play any significant role under our experimental condi- 
tions. The exponential solution to Eq. [2] for t > is 
N(t) = JVoo(l - e~ 7wt ), where y tot = y Na + y bg is the to- 
tal trap loss rate. Each loading curve was fit to this func- 
tional form to extract Noo and y tot . The initial loading rate 
is then determined as the product of the two extracted pa- 
rameters R = Naajtot- 

As elucidated by [22], for a vapor cell MOT both the cap- 
ture rate R and the loss rate are proportional to the 
density of sodium and therefore to each other. Thus we 
can parameterize R = ay^ a for some a. From this the 
extracted fitting parameters take on specific functional 
forms 



R 

7tot = - + 7b g 
a 



(3) 



and 



R aR 

Nee = — = ~ (4) 

7tot R + oTbg 

In Figure [2] we test the validity of the above equations 
by plotting the experimentally determined value of y to t 
and Nee versus R, which was varied by changing the in- 
tensity of the LIAD using the blue LED. In Fig|2ja), the 
data clearly shows a linear relationship between y tot and 
R whose slope yields a = 4.4 x 10 7 , and whose intercept 
yields ybg = 0.19 s . Moreover, the steady state atom 
number fits well to Eqn. (0 with the same value of a and 
7b g within statistical uncertainty. 

The parameter a represents an upper bound for the atom 
number in the MOT. In the limit where sodium collisions 
dominate the MOT loss rates, — > R/y^a = oc becomes 
independent of vapor pressure and depends only on the 
ratio of thermal to capture velocities for the atom, the vol- 
ume enclosed by the trapping beams and upon other opti- 
cal parameters (laser intensity, detunings and repumping 
sideband strength). In this situation, further increasing R 
by turning up the light intensity used for LIAD will only 
shorten the time taken to reach steady state, not increas- 
ing the maximum trapped atom number. In Figure Q] one 
can clearly see that for high LIAD intensity TVoo, which 
appears to be saturating near 4 x 10 7 . The inset shows, 
however, that the loading rate R continues to increase lin- 
early with LIAD intensity up to at least 4mW/cm 2 . The 




FIG. 2: Atom number and loading rate dependencies. By fitting 
each curve in Fig. [T] to the solution of Eq. 10, the total loss 
rate y tot (a), and the total number of captured atoms N m (b) are 
determined and plotted against the loading rate R = ytotW»- The 
solid lines are fits from the Eqns. l(3j and 10. 



maximum measured R was almost 50 times larger than in 
the absence of LIAD, while the steady state atom number 
was a little less than 10 times larger. 



A. The Role of Background Gas 
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FIG. 3: Influence of background gas collisions. We present 
three data sets of the normalized atom number v for two differ- 
ent vacuum regimes. The two UHV data sets were obtained by 
varying the blue LED voltage after different elapsed times since 
the Na dispenser was fired, while the HV data set was obtained 
by changing the dispenser current. The solid curves result from 
the direct use of Eqn. l[5j, with ybg determined by Eqn. (O. 

To further confirm the model of 12211 and clarify the role 
of background gas collisions, it is useful to normalize the 
steady state atom number to its maximum value a. If 
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FIG. 4: LIAD loaded MOT is compatible with long lifetimes. 
The entire loading and decaying dynamics of the MOT with the 
blue LED LIAD is shown. The LED was turned on at time 
t = 65 s, and turned off at t = 105 s. By fitting the remaining 
trapped atom number for / > 105 s we extracted the loss rate 
of 0.10 s" 1 , which is comparable to the value inferred from the 
steady state properties of the model Eqn. {2j, as well as with the 
background MOT loaded previously, between 5 and 65 s. 

we define a dimensionless variable v = Nco/a it can be 
expressed via Eqn. (0]) as 



7Na + 7bg 

TNa is proportional to the density of sodium which can be 
varied by changing the LIAD intensity or the dispenser 
current. Therefore, v runs from 0, when no sodium 
is present (in the background vapor), to 1, when col- 
lisions with Na completely dominate the MOT losses. 
The crossover between these two regimes occurs when 
TNa = Tbg, the background gas collision rate. We changed 
the vacuum regime in order to study the dependence of 
this parameter on the MOT loading dynamics. 

For both high vacuum (HV) and ultra-high vacuum 
(UHV) conditions, loading curves were acquired and an- 
alyzed in a manner similar to Figure Q] Thus a was com- 
puted and both v and y^a = R/a could be determined 
experimentally. The result is plotted in Figure [3] We 
observed a dramatic dependence of the measured quanti- 
ties on the vacuum conditions. The data sets all obey the 
same functional form given by Eqn. ((5]), but the crossover 
point occurs for much lower y^ a , i.e., sodium vapor den- 
sity, for UHV when compared to HV. 

From a fit to Eqn. ([3]) we can extract the values of 
the background gas collision rate. Under HV condi- 
tions where the pressure was 6 x 10~ 9 torr, we obtained 



7bg ~ 1.8 s -1 . Under UHV conditions (pressure be- 
low 1 x 10" 11 torr), the background gas collision rate 
was found to vary depending on the time elapsed since 
the last firing of the dispenser. The UHV curve with 
7bg = 0.19 s _1 was taken a few hours after the firing 
whereas the curve with y^,„ = 0.13 s _1 was taken a few 
months after the last firing. The ratio (ybg)Hvl{jbg)uHV ~ 
10. By contrast, the vacuum gauge measured a pressure 
ratio of about 600. This is not an uncommon feature of 
UHV systems: due to mounting constraints the vacuum 
gauge had to be installed closer to the pumps where the 
pressure is typically lower than at any other location in 
the vacuum chamber, including the atom trap. In fact, the 
data in Figure [2ja) is more reliable to measure the life- 
time of the trapped atoms than the pressure gauge read- 
ings. 

The MOT saturates at fairly low sodium vapor density in 
the UHV regime, and it is typically larger and brighter. 
Under UHV conditions beyond the crossover point, in- 
creasing the vapor density by a factor of 5 increased 
the MOT atom number by only 14%. This observation 
demonstrates that very low Na background vapor density 
is needed to completely fill the MOT at an appreciable 
atom number of 3.7 x 10 7 . Moreover, the UHV MOT did 
not require the dispenser to be fired over the course of 
many months, when running it no more than three to four 
times a week fl4£] . 

The collisional loss rates in the range y tol — 0. 1 to 0.2 s _1 , 
inferred from the data in Figj3] suggest that the LIAD 
loaded MOT should be compatible with evaporative cool- 
ing to BEC. We have also measured the time constant 
directly by observing the MOT dynamics. The result is 
plotted in FigHl where we recorded the time evolution of 
the atom number before, during and after the LED was 
switched off. In that data set, the MOT atom number de- 
cayed from a peak of 10 7 around t — 105 s to an asymp- 
totic value of about 2.5 x 10 6 atoms with a 1/e time con- 
stant of t = 10 s. 

In the absence of LIAD, a small background Na partial 
pressure exists, allowing for a small MOT to be filled be- 
fore t = 65 s. When the LED is turned on it creates a tran- 
sient Na vapor pressure of much higher density that is re- 
sponsible for the spike in the trapped atom number. Upon 
shutting off the LED the extra vapor becomes quickly ex- 
tinguished by being adsorbed onto the cell walls. Colli- 
sions with this background sodium vapor as well as resid- 
ual (non-sodium) vapor are the cause of the observed 10 
second decay time. 
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B. Mechanism of desorption 

According to Yakshinskiiy and Madey fT3l . the photo- 
desorption of neutral sodium is a charge transfer process 
for Na + ions bound to the silica (SiCh) surface. An elec- 
tron can be excited from SiC>2 by photon bombardment. 
Above a threshold energy this electron transfers to the un- 
filled Na 3s level. The neutral Na atom then finds itself on 
a repulsive potential causing it to quickly dissociate from 
the surface. Previous work demonstrated a steep depen- 
dence of the neutral sodium yield on photon energy in the 
range of a few eV using a mercury arc lamp. However, 
the large error bars caused by the broad wavelength spec- 
trum of the mercury lamp precluded a more precise de- 
termination of the wavelength dependence in the region 
of the threshold. Fig.l of reference [12] shows an error 
bar spanning the range of approximately 2.0 - 3.8 eV for 
the lowest energy filter setting. No additional informa- 
tion is available within or below that energy range as the 
experiment could not detect any appreciable signal. For a 
more precise determination of the desorption due to solar 
radiation incident on planetary bodies such as the Moon, 
it is the photon flux from the visible and ultraviolet radi- 
ation that plays the most important role. Therefore there 
is a clear need to examine this region in greater detail. 

To address this issue we recorded the MOT loading 
curves using eight independent LEDs whose wavelengths 
ranged from 590 nm to 365 nm. A black paper mask was 
custom built to fit over a small part of the Pyrex vacuum 
cell with a 2 by 2 cm aperture on both front and back sur- 
faces. This approach allowed us to precisely define and 
control the area exposed to the LIAD photons. Before 
the runs, the LEDs were individually adjusted to provide 
2 mW/cm 2 of light irradiance through the mask aperture. 
This was experimentally found to optimize the LIAD sig- 
nal to noise ratio over the range of wavelengths used. 
Since the spatial properties of each LED were slightly 
different, we carefully measured their individual intensity 
profile within the exposed area and applied an area aver- 
age correction factor to the data. The photon energies and 
energy uncertainties are provided by the manufacturer's 
specifications sheets 11231 12411 . 

To acquire the data, each run consisted of two full MOT 
loading curves, similar to those in FigQ] taken in two 
steps. In the first step, the LED was kept off and the MOT 
was fully loaded from the background sodium vapor for a 
period of about 40 s. The MOT current was then immedi- 
ately turned off and the LED was turned on and allowed 
to equilibrate thermally for about 20 s. In the next step the 
MOT was reloaded with the LED on for 40 s. Thus, in the 
first step, the MOT was loaded from the background Na 
vapor alone, while in the second step both background 



vapor and photo-desorbed atoms were present. The final 
result was then computed as the difference between the 
rates dN/dt\t=a measured with LEDs on and off, yielding 
the net loading rate due to photo-desorption. 
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FIG. 5: Photon energy dependence of desorption of sodium 
atoms from Pyrex glass. Plotted is the excess MOT loading rate 
due to LIAD for an average intensity of 2 mW/cm 2 . Each point 
is the average of 3-4 MOT loading curves. The inset shows the 
same data plotted on a logarithmic scale with the y-error bars 
resulting from statistical data analysis. The energy error bars 
correspond to the typical half-width uncertainties directly pro- 
vided by the manufacturer's specifications sheets. 

The results are shown in Figure|5] where one can clearly 
see that the desorption rates become significant only for 
photon energies above 2.6 eV (470 nm). For longer wave- 
lengths we are unable to detect any significant LIAD ef- 
fect, consistent with a threshold phenomenon. In the ab- 
sence of a theoretical model we have refrained from fit- 
ting our data to any functional form. Nevertheless, we 
are able to detect a significant sodium desorption in the 
range 2.6 to 3.4 eV, a region below the detection thresh- 
old of previous work [12]. Part of the improvement lies 
in the use of LEDs rather than a Hg lamp, which has nar- 
rowed the spectral bandwidth associated with each data 
point by an order of magnitude on average. More sig- 
nificantly, however, our atom trapping technique is very 
sensitive to small numbers of desorbed atoms. For exam- 
ple, the majority of the data we presented in Figures [T]|4] 
were taken using light at 455 nm (2.72 eV), very close to 
the threshold observed for the LIAD effect on sodium. 
Therefore, our data on Pyrex complement the studies on 
pure Si02 surfaces lfl2ll25ll . and there could be variations 
in the material surface properties that may affect the PSD 
as well. 

The MOT loading rate depends upon both the photo- 
desorption cross-section, as well as the MOT capture ef- 
ficiency. If the desorbed atoms re-thermalize with the 
cell walls before they enter the MOT laser beams, as we 
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believe to be the case for our experimental scenario, the 
latter should be independent of the wavelength of the des- 
orbing light. In this limit the loading rate of Figure |5]be- 
comes a relative measure of the photo-desorption cross- 
section. However, it is important to note that the photo- 
desorption process does not necessarily need to produce 
a Maxwell-Boltzmann distribution of velocities, and the 
MOT capture efficiency need not be wavelength indepen- 
dent, as is often assumed to be the case. 

Larger photon energies would increase the desorption 
rate dramatically, an advantage when producing large 
MOTs very quickly; however, shorter wavelengths will 
be significantly attenuated through the Pyrex walls. Fu- 
ture work could explore the difference between glass and 
quartz, although it is questionable whether the latter will 
be as useful for producing larger MOTs [8]. We also plan 
to examine porous materials more closely resembling the 
lunar surface. 



IV. CONCLUSION 

We have demonstrated that LIAD can be used effectively 
for fast loading of a sodium vapor-cell based magneto- 
optical trap, allowing for rapid control over the alkali 
background vapor pressure inside the cell. The same 
desorption process plays an important role in planetary 
science, albeit at drastically different temperatures from 
a laboratory BEC experiment. Studies by the surface 
science community have increasingly beg un to quantify 
radiation-induced chemical changes Jl2l |25L 12611 . The 
high sensitivity (potentially single atom) of magneto- 
optical traps may be a powerful tool for surface science. 

C.R. thank Thomas Orlando for useful discussions. This 
work was supported by the U.S. Department of En- 
ergy, the Army Research Office and the National Science 
Foundation. 
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